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Abstract 

A vertex model introduced by M. Bowick, P. Di Francesco, O. Golinelli, and E. 
Guitter [Nucl. Phys. B 450, 463 (1995)] describing the folding of the triangular lattice 
onto the face centered cubic lattice has been studied in the hexagon approximation 
of the cluster variation method. The model describes the behaviour of a polymerized 
membrane in a discrete three-dimensional space. We have introduced a curvature 
energy and a symmetry breaking field and studied the phase diagram of the resulting 
model. By varying the curvature energy parameter, a first-order transition has been 
found between a flat and a folded phase for any value of the symmetry breaking field. 
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The macroscopic behaviour of membranes fluctuating in the euchdean space crucially 
depends on the microscopic characteristics of the system [jl| . For example, a fluid membrane 
without self-avoiding interactions is always expected to be in a crumpled state 0, indepen- 
dently from the stiffness attributed to the membrane, while a rigid phantom polymerized 
membrane, which is a network of molecules with flxed connectivity, is expected to be stable 
in a flat phase The existence of different classes of membrane systems suggests to con- 
sider speciflc lattice models paradigmatic for each class, which can be useful for analytical 
and numerical calculations. While many-component membrane systems can be sometimes 
expressed as Ising-like models [^, it is more difficult to describe the statistical behaviour of 
a single membrane in terms of usual lattice models with a local hamiltonian. In Ref . , a 
D = 2 vertex model has been introduced which describes the behaviour of a single phantom 
polymerized membrane with bonds of flxed length embedded in a discrete d = 3 space. The 
aim of this communication is to study the phase diagram of the membrane model of j^. 

Models of polymerized membranes were introduced in |^, and studied using a Landau- 
Ginzburg evaluation of the energy of the system in |^, |rU|. From the above studies, 
and from numerical simulations []TT|, it comes out that, when excluded volume effects are 



not considered, by varying the strength of a bending energy term which favours flat 
conflgurations, a critical transition arises separating a flat phase at large K from a crumpled 
phase at small K. 

Generally, the length of the bonds of a polymerized membrane can vary being subjected 



to elastic forces H]. In [12| the simple case of a triangular network with bonds of flxed length 



embedded in a (i-dimensional space has been flrst considered. This folding problem has been 
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studied in |T^, |l^, ^ in the case of a two-dimensional embedding space; here the normals 
to the triangles of the network can point only "up" or "down" in some direction. This model 
can be mapped on a 11-vertex model equivalent to a constrained Ising model with some spin 



configurations forbidden |13]; a first-order transition has been found to occur between a fiat 



and a disordered phase [14, 15 . 



The more complicated problem of the folding of the triangular lattice in a three-dimensional 
embedding space has been formulated in with the embedding space discretized as a face 
centered cubic lattice. In this model the plaquettes of the triangular lattice are mapped, by 
folding, onto those of a face centered cubic lattice, so that two adjacent plaquettes form an 
angle which can take up only four different values (see Fig. 1). In the following we introduce 
a term representing the stiffness of the network and study the phase diagram of the resulting 
model by applying the cluster variation method (CVM) in the same approximation used in 



15 . 



The vertex model of can be defined as follows. Two Z2 variables, named Zi and (Xj, 
are assigned to each plaquette of the triangular lattice, and hence to each site i of the dual 
hexagonal lattice. The relative values of these variables on adjacent plaquettes (say 1 and 2) 
determine the angle that is formed by the plaquettes, according to the following rules: for 
zi = Z2 and o"! = 0-2 we have no fold between the plaquettes; for Z2 = —zi and a2 = cri we 
have an acute fold, with an angle 6 such that cos 6 = 1/3; for Z2 = —zi and a2 = — ci we 
have an obtuse fold, with cos 6' = —1/3; finally, for Z2 = z\ and 02 = — o"i we have a complete 
fold, with the plaquettes lying one on top of the other (see Fig. 1). 

It has been shown in that the variables Zi and cTj have to satisfy two constraints, or 



folding rules, in order to describe a folding configuration over the face centered cubic lattice. 
Such folding rules take the form 

6 

^ CTj = mod 3, (1) 

i=l 

where the index i runs around a hexagon, and 

E^^l^ A.c = 0mod2 c=l,2 , (2) 
i=i ^ 

where z-j = Z\ and 

A.e=|i ifp^i^.=^-od3 ,^1^^ 6 . (3) 

' I otherwise ' ' ' ' ^ ' 

We shall assume that a fold between two adjacent plaquettes has an energy cost, due 
to curvature, given by —Kcos6, where 6 is the angle between the normal vectors to the 
plaquettes. In terms of our Ising-like variables such a term can be rewritten as —Kai(jj{l + 
2ziZj)/3, manifestly symmetric under the global transformations Zi — > —Zi,Wi and Uj 
—ai,'^i. We will introduce also a term which breaks this symmetries, the analog of the 



magnetic field in the ordinary Ising model. In the two-dimensional case |T^, |I5| it was 
quite easy to define a symmetry-breaking field coupled to the direction of the normal to a 
plaquette, since there the spin variable associated to each plaquette denoted whether the 
normal was pointing up or down. In the present case and with the present parametrization 
this is no more possible, since the pair of values {zi,ai} does not determine the orientation 
of the plaquette. However, we are not interested in a true magnetic field, which cannot be 
given a physical meaning in this model, but just in a symmetry-breaking term, which can 
be realized in several ways. A simple choice, which favors only one of the four possible fiat 
states (a fiat state is easily seen to be characterized by the condition {zi, a^} independent of 



i), leaving on the same ground the remaining three, is —hSz^, 16^^,1- There are of course other 
possibihties, and in the following we shall consider one of these in order to show that the 
phase diagram is, roughly speaking, qualitatively independent of this choice. 

We are thus led to consider the following hamiltonian (energies are given in units of k-^T), 
defined on the dual hexagonal lattice: 

n = - + 2^*^i) - (4) 

where the first sum is over nearest neighbor pairs. 

As for the two-dimensional model, the phase diagram will be investigated by means of the 
hexagon approximation of the cluster variation method (CVM), which has been thoroughly 



described in ||T5[. This requires to introduce a hexagon density matrix p6({-2i, cTj}), indexed 

by the hexagon configurations, a link density matrix p2 and two site density matrices piA and 

PiBi one for each of the two interpenetrating sublattices of the hexagonal lattice. Because 

of the folding rules Eqs. |1| and ^ one has not to consider 2^^ = 4096 hexagon configurations, 

but only 384, and hence 384 pg elements (this number may be slightly reduced by taking 

into account degeneracies, but this would lead to a more involved formulation). The link 

and site density matrices can be defined as partial traces of pe and, in analogy with [|T^], we 

can write the variational free energy density (in units of k-oT) 

K h 
f = -—TT[{aia2{l + 2ziZ2)p2]--[piA{z = l,a=l) + piB{z = l,cr=l)] 

13 1 
+ 2Tr(P6lnp6) - -Tr(p2lnp2) + - [Tr(piA Inpi^) + Tr(piB Inpi^)] (5) 

+A(Tr p6 - 1), 

where A is a Lagrange multiplier which ensures the normalization of pe (and hence also of p2 



and PiaiPib)- This variational functional must, in general, be minimized numerically, and 
this can be done easily by an iterative method, as explained in [p!5[] . 

At infinite temperature (or vanishing K and h) we obtain the entropy (per site) = In g, 
where q ~ 1.42805, in very good agreement with the transfer matrix estimate q = 1.43(1) 
found in 0, although slightly lower than the best lower bound 1.43518 obtained from the 
analysis of two-dimensional folding in a staggered field 0. Furthermore we obtain {zi)A = 
{zi)B = and ((Tj)^ = —{o'i)B — 0.87456, which, together with the values of the link 
density matrix elements, indicate a marked preference of the triangular lattice for obtuse 
and complete folds. 

Let us now describe the main features of the phase diagram of the model Eq. H, which 
has been obtained by means of the CVM and reported in Fig. 2, in the K — h plane. The 
solid line separates the fiat phase (high values of K and /i, or low temperature) from the 
disordered, or folded, one (low values of K and h, or high temperature). The transition is of 
first order on the whole line; on the = axis the transition occurs at i^' = 0.18548, while 
on the K = Q axis it occurs aX h = 0.49839. The two phases can be distinguished by means 
of the energy-like correlation function (crjcrj(l + 2ziZj)/?)), which is negative in the folded 
phase and saturates to 1 in the fiat phase, indicating that the plaquettes are all parallel. 
As a consequence, in the fiat phase the entropy vanishes and the free energy reduces to the 
internal energy contribution, 

3 

h^t = --K-h. (6) 

The phase diagram of Fig. 2 is then easily obtained by comparing the value of the free 
energy of the disordered phase with /flat above. 

6 



The first point to be discussed is that enlarging the embedding space has not turned the 
phase transition into a second order one, in contrast with what one could expect 0. On 
the contrary, the first order character of the h = transition seems to be enhanced in the 
three-dimensional case, as suggested by the jump in the energy-like correlation function, 
which is about 1.237 in the present case against 1.047 in the two-dimensional case. These 
features should be described correctly by the CVM approximation, as suggested by the good 
agreement with the transfer matrix results obtained both in the two-dimensional case W5 
and, for the entropy at infinite temperature, in the present case. 

The curvature term in Eq. ^ without the constraints Eqs. corresponds to the 



anisotropic Ashkin-Teller model with a particular choice of the parameters. We have 
studied in the CVM pair approximation the Ashkin-Teller model with the hamiltonian 
'Hat = ~K/3J2<ij> o'i'^ji^ + '^ZiZj) and we have found two critical transitions at Kd = 0.824 
(where the symmetries z ^ —z and a —a are separately broken, but their product is 
preserved) and Kc2 = 1.648 (where even the product symmetry gets broken). In analogy 
with the result obtained in the two-dimensional case [T^, we see that the introduction of 
defects which relax the constraints is necessary to smooth the transition. 

Finally, it has already been observed before that the introduction of a symmetry-breaking 



field in the hamiltonian is not a trivial task as in the two-dimensional case [T5[ . It is therefore 
worth asking whether the basic features of our phase diagram depend on the choice of this 
field. For this reason we have considered a second model by introducing the /i-field term in 
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the hamiltonian in the following way 

Kit = E + 2'2i'2j) - ^ E ^^(1 + 2-^0, (7) 

where h is non-negative. This amounts to let all the plaquettes interact, with an energy h, 
with a "ghost" plaquette which is fixed in the state characterized by (z = 1, cr = 1). In Fig. 
2 it is shown the phase diagram of the model Eq. |^ as well; the transition line is the broken 
one and is still a first order one. It is clear that the main features of the phase diagram have 
not been changed by replacing the old /i-field term by the new one. 

To summarize, we have studied in the CVM hexagon approximation the phase diagram 
of a vertex model which describes the folding of a triangular network onto a fee lattice, 
subjected to a bending energy and a symmetry-breaking field. We have determined the 
folding entropy at infinite temperature, which is in very good agreement with a previous 
transfer matrix estimate, and the phase diagram of the model, which turns out to be quali- 
tatively independent on the choice of the symmetry-breaking field. The fiat and the folded 
phases are separated by a first-order transition, which is even stronger than that obtained 
with a two-dimensional embedding space. The lack of critical behaviour, which is expected 
for a polymerized phantom membrane, may be due to the particular choice of the discrete 
embedding space. 

One of us (G.G.) thanks Amos Maritan for a discussion about the subject of this paper. 
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Figure Captions 



Figure 1: The four possible foldings of two adjacent plaquettes of the triangular lattice 
embedded in a fee lattice are shown. From the left to the right and from the top to the 
bottom: no fold, acute fold (70°32'), obtuse fold (109°28') and complete fold. The dots 
represent the vertices in an elementary cell of the fee lattice. 

Figure 2: Phase diagram in the K — h plane: the sohd and dashed hues correspond to 
hamiltonians H and Haxt, respectively. 
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